Orbital-selective magnetism in FeAs-based superconductors 
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We report first-principles studies of lanthanide-series (Ln) iron oxypnictides LnFeAsO (Ln=La, 
Ce, Pr, Nd, Sm, and Gd) which show two competing orbital-selective magnetic phases in small-Fe- 
moment regime: an in-plane d xy magnetic phase, itinerantly driven by orbital selection of Fermi- 
surface nesting, and an out-of-plane d yz magnetic phase, driven by local interactions. The Fe 
magnetic moments in the two phases show different coupling strengths to Fermi-surface electrons 
orbital-selectively, suggesting different roles in superconductivity and in antiferromagnetism, and 
making orbital characters of the moment resolvable by measuring the electronic structures. These 
results show that orbital sensitivity is a key factor in the magnetic properties of FeAs-based materials. 

PACS numbers: 71. 15. Mb, 71.20.-b, 74.70.Xa, 75.25.-j 



LaFeAsOi-^Fz [l[ and related compounds show un- 
conventional superconductivity_(SC) in the vicinity of 
antiferromagnetism (AFM) [2h16| . Among various fami- 
lies of iron pnictides and chalcogenides, lanthanide-series 
(Ln) iron oxypnictides LnFeAsO reach the highest super- 
conducting transition temperature (T c ) with doping. Re- 
ported T c in doped LnFeAsO increases dramatically from 
26 K up to 55 K with Nd or Sm substitution for La and 
then T c decreases slightly in doped GdFeAsO [H 0, H, [ll|- 
15j . With this strong variation of T c , LnFeAsO is suited 
for studying material dependence of T c . Since FeAs- 
based materials are featured with multiple Fermi surfaces 
(FSs) with strong orbital characters, many studies have 
been done on FS nesting, local-moment interactions, and 
orbital orderings [l7-27|. However, no specific feature 
has been pinpointed in electronic and magnetic proper- 
ties correlated with material dependence of T c except for 
As positions and their effect on FS. 

In this Letter, we report, for the first time, orbital- 
selective magnetic phases in FeAs-based materials found 
by first-principles density functional theory (DFT) cal- 
culations. Two phases are identified with distinctive or- 
bital characters of Fe magnetic moments, d xy and d yz , 
which originate from FS nesting and local interactions, 
respectively. The two phases, while competing with each 
other, show different electronic band structures caused by 
orbital-selective coupling strengths between Fe magnetic 
moments and FS electrons. Our results suggest that spin 
or magnetic interactions for unconventional SC should 
have strong orbital sensitivity in FeAs-based materials. 

Our calculations are based on ab initio norm- 
conserving pseudopotentials and the Perdew-Burke- 
Ernzerhof-type generalized gradient approximation, as 
implemented in the SIESTA code [281 ] . During self- 
consistent iterations, constraints are imposed on the elec- 
tron spin density to reduce the size of the Fe magnetic 
moment (tofc) to specific values smaller than uncon- 
strained results. We assume an oxidation state of +3 
for all considered Ln elements, treating their 4/ orbitals 
as core orbitals 111, 3(|. Experimental lattice constants 
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FIG. 1: (Color online). Orbital characters of Fe magnetic 
moments: (a) LaFeAsO and (b) SmFeAsO when rriFe is 0.15 
/xb per Fe atom. Each atomic structure is a y2 x \/2 x 1 
supercell with four Fe atoms in the single-stripe-type AFM. Fe 
magnetic moments, in dark and light grays (purple and light 
gray), are ordered antiferromagnetically along the x axis and 
ferromagnetically along the y axis. Calculated Fe magnetic 
moments have two distinct orbital shapes, d yz in LaFeAsO 
and d xy in SmFeAsO. The considered mp e of 0.15 /j,b is about 
one half of experimental values, 0.36 fiB in LaFeAsO [j| and 
0.34 fj, B in SmFeAsO [3. 



and atomic positions at the high-temperature tetragonal 
phase are used for the calculations 0, 0, EH lil 31|. 

An intriguing new finding in our present work is Ln- 
dependent distinctive d-orbital characters of the Fe mag- 
netic moment in LnFeAsO in small mp c regime. Fig- 
ure 1 shows two prototypical cases of LaFeAsO and Sm- 
FeAsO. For mp c of 0.15 Bohr magneton (/is) per Fe 
atom, LaFeAsO shows out-of-plane d yz character, but 
SmFeAsO shows in-plane d xy character. In previous 
studies, the main concern was not d xy versus d yz , but 
d yz versus d zx regarding symmetry lowering from tetrag- 



onal to orthorhombic structures. In our work, symmetry 
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FIG. 2: (Color online). Orbital analysis of Fe magnetic mo- 
ments in LnFeAsO. (a)-(f) Orbital contribution to the Fe mag- 
netic moment as a function of rriYe per Fe atom in LnFeAsO. 
Downward arrows indicate critical values of rriYe at which or- 
bital characters switch between d xy and d yz . In LaFeAsO, 
the critical mr e is very close to zero, (g) Critical m^e versus 
the As height from the Fe layer, (h) Ln-dependences of criti- 
cal m.Fe and experimental T c . Experimental T c 's are 26 K in 
LaFeAsO0.s9F0.11 Q, 41 K in CeFeAsO0.84F0.i6ji, 52 K in 
PrFeAsO0.89F0.11JP], 54.3 K in NdFeAsOi-^ GjL 55 K in 
SmFeAsOo.gFo.i BjL and 54 K in GdFeAsOi_„ fijl. 
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FIG. 3: (Color online). Orbital characters and pocket areas 
of Fermi surfaces (FSs) in nonmagnetic LnFeAsO. (a) FS of 
nonmagnetic SmFeAsO at k z = and (b) at k z — ir/c. The ki 
and &2 axes are along reciprocal lattice vectors of the nonmag- 
netic unit cell having two Fe atoms. In (a), the 45°-rotated 
inner square is the first Brillouin zone of the y/2 x y/2 x 1 su- 
percell having four Fe atoms, where k x and k y axes are along 
the x and y axes in Fig. 1(a), respectively. Three hole pock- 
ets (a, ft, and 7) are located at the center and two electron 
pockets (5 and e) at the corner. Orbital characters, d xy , d yz , 
and d zx , are represented in blue, red, and green, respectively. 
The a, 7, and S pockets are derived from d yz and d zx , the j3 
pocket is from d xy , and the e pocket is from d xy , d yz , and d zx . 
(c) FS pockets at k z — 0. The 5 and e pockets are plotted as 
closed loops to highlight changes in orbital characters around 
the pockets, (d) Areas of the five FS pockets at k z = and 
(e) at k z — tt/c as a function of Ln elements. 



is lowered by imposing the single-stripe-type AFM, yet 
d xy character is predominant in SmFeAsO. 

For systematic study, we analyze orbital contributions 
to the Fe magnetic moment in LnFeAsO as a function of 
rriFe [Figs. 2(a)-(f)]. We find that either d xy or d yz char- 
acter is dominant, switching abruptly at a critical value 
of rriFe in each compound. The critical rriFe is 0.014 /ib in 
LaFeAsO and increases to 0.23 /xb in GdFeAsO. Except 
for GdFeAsO, the d xy (d yz ) character appears at smaller 
(larger) mp c than the critical value. In LaFeAsO, the 
critical tof c is so small that the d yz character is domi- 
nant [Fig. 2(a)]. In GdFeAsO, unlike other compounds, 
the d yz character is larger than the d xy character at m-pe 
smaller than the critical moment [Fig. 2(f)]. 

The overall increase of the critical rriFe from LaFeAsO 



to GdFeAsO is correlated most strongly with the increase 
of the As height from the Fe plane among various struc- 
tural parameters [Fig. 2(g)], and correlated with experi- 
mentally reported T c [Fig. 2(h)]. This suggests that the 
d xy magnetic moment is closely related with measured 
T c in the doped LnFeAsO for Ln from La to Sm, while 
GdFeAsO shows different and more complicated orbital 
features possibly because GdFeAsO is beyond an optimal 
chemical composition for T c [J 0, @, [HHS] • 

To understand the origin of the orbital-distinctive mag- 
netic moments, we analyzed orbital characters of FS in 
nonmagnetic LnFeAsO. As a prototype, SmFeAsO has 
three cylindrical hole pockets [a, (3, and 7 in Figs. 3(a)- 
(c)] and two electron pockets [S and e in Figs. 3(a)-(c)] 
at the Fermi energy (Ep). Orbital analysis shows that 
the a, 7, and S pockets are derived from Fe d yz and d zx 
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tions of d yz and d zx 



moment of d yz or d 



orbitals, the f3 pocket is from d xy , and the e pocket is 
from d xy , d yz , and d zx [Figs. 3(a)-(c)]. For FS nesting, 
we consider six pairs of hole and electron pockets: a-S, 
j3-5, 7-5, a-e, (3-e, and 7-e. Among these, a-S nesting is 
not significant because the two pockets have quite dif- 
ferent pocket areas in all LnFeAsO [Figs. 3(d) and (e)], 
and (3-5 nesting is not effective because orbital difference 
of the two pockets degrades effects of nesting greatly. In 
addition, j-S and a-e nesting effects should be very weak 
because of orbital mismatch due to out-of-phase alterna- 
characters around the pockets. Thus, 
we need to examine only /3-e and 7-e nesting regarding 
the orbital character of the Fe magnetic moment. 

Since the /3-hole pocket is derived from the d xy orbital, 
[3-e nesting may induce Fe magnetic moment of d xy char- 
acter. Similarly, 7-e nesting may induce Fe magnetic 
character. Calculated FS pocket 
areas as a function of Ln elements [Figs. 3(d) and (e)] 
show that the /3-hole pocket area is sensitive to Ln ele- 
ments, approaching the e-electron pocket area for heavier 
Ln elements, while the 7-hole pocket area is almost insen- 
sitive to Ln elements and close to the e-electron pocket 
area. This indicates that (3-e FS nesting is a key factor in 
the electronic structures in LnFeAsO which changes the 
orbital character of the Fe magnetic moment. 

Although it is not involved in the FS nesting, the area 
of the 5-electron (d yz and d zx ) pocket also depends sen- 
sitively on Ln elements [Figs. 3(d) and (e)]. With heav- 
ier Ln elements, the /3-hole (d xy ) band moves upward 
in energy and the 5-electron band moves downward in 
energy, resulting in electron transfer from the d xy or- 
bital to the dyz and d zx orbitals. Since the 5-pocket is 
not related with FS nesting, the increased d yz and d zx 
electrons can contribute to magnetic moments only by 
non-FS-nesting mechanism, e.g., local- moment interac- 
tions 0, [3]. Therefore, while the enhanced d xy mag- 
netic moment from LaFeAsO to SmFeAsO is due to bet- 
ter /3-e FS nesting, the enhanced d yz magnetic moment 
in GdFeAsO is due to increased d yz electrons in the S- 
pocket which favor non-FS-nesting mechanism. 

To find out effects of different orbital characters of 
Fe magnetic moments on the electronic structures, we 
obtained orbital-resolved band structures and projected 
density of states (PDOS) in SmFeAsO, as a prototype of 
LnFeAsO, in nonmagnetic and magnetic phases (Fig. 4). 
With zero m,Y e , energy bands have distinct orbital char- 
acters along high-symmetry lines, yielding PDOS slowly 
varying near Ep [Figs. 4(a) and (b)]. 

With the d xy magnetic moment, strong anti-crossing 
occurs between d xy bands, opening orbital-dependent en- 
ergy gaps at Ep [Fig. 4(c)]. In the TX line, d xy and d zx 
bands repel each other at Ep and at -260 mcV. In the TY 
line, d xy bands open a large energy gap near Ep [marked 
with A xy in Fig. 4(c)], while d zx bands are intact. Except 
for anti-crossing at -100 meV in the TY line, d yz bands 
are almost unchanged. The d xy PDOS is greatly modi- 
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FIG. 4: (Color online). Orbital-resolved electronic band 
structures and projected density of state (PDOS) in Sm- 
FeAsO: (a) and (b) for nonmagnetic phase, (c) and (d) for the 
d xy phase (mF C = 0.15 (j,b), and (e) and (f) for the d yz phase 
(mF C = 0.20 /ib). Band structures are plotted along high- 
symmetry lines in the first Brillouin zone of the y/2 x V2 x 1 
supercell for the single-stripe-type AFM, where M = (0, tt), 
X = (it/2, it/2), T = (0, 0), and Y = (-tt/2, tt/2) with respect 
to the unit-cell reciprocal axes, k\ and k-z, shown in Fig. 3(a). 
Blue, red, and green represent d xy , d yz , and d zx orbital char- 
acters, respectively, and Ef is set to zero. Horizontal arrows 
indicate extreme band energies at Y. In (c), A xy indicates 
energy gap between d xy bands in the TY line. In (e), A yz 
indicates energy gap between d yz bands in the FY line. 



fied with a partial-gap opening at Ep and a huge peak 
right above E F [Fig. 4(d)]. The d zx PDOS also shows a 
reduction near Ep because of coupling to d xy [Fig. 4(d)]. 
Despite these large effects of the d xy magnetic moment on 
the electronic structures near Ep , deformation of bands 
is mostly near the crossing points, with no shift of origi- 
nal band edges [Fig. 4(c)]. This confirms that the origin 
of the d xy magnetic moment is due to FS nesting. 

In contrast, with the d yz magnetic moment, d yz bands 
are deformed in the whole Brillouin zone much more 
greatly than anti-crossing of bands [Fig. 4(e)]. Whole 
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upper part of d yz bands and whole lower part of d yz 
bands are expelled from each other [Fig. 4(e)], so even 
the top of the d yz hole band is pushed up by 100 meV 
and the bottom of the d yz electron band is pushed down 
by 100 meV. This change of d yz bands in the whole Bril- 
louin zone indicates that the driving mechanism for the 
d yz magnetic moment is not FS nesting (local in k-space) 
but lowering of the total energy of electrons by formation 
of local magnetic moments (local in real space). Despite 
the large change in the d yz bands, all the other bands are 
almost unchanged except for some anti-crossings. 

Since different orbital characters of the Fe magnetic 
moment result in orbital-dependent difference in elec- 
tronic structures, the orbital characters of the Fe mag- 
netic moment can be distinguished by, e.g., angle- 
resolved photoemission spectroscopy of detwinned sam- 



ples [27J, 1221 ■ A s established by previous studies, DFT 



electronic band structures are qualitatively consistent 
with experimental results in FeAs-based materials al- 
though the DFT results overestimate experimental band- 
widths by a factor of 2 to 3 near Ep. Thus, detailed com- 
parison of our band structures with experimental results 
will reveal orbital characters of the Fe magnetic moment, 
and thereby their roles in AFM. 

As shown above, the d xy magnetic moment is driven by 
orbital-selective FS nesting, while the d yz magnetic mo- 
ment is related with spatially localized interactions such 
as local-moment interactions. Thus, the competition of 
d X y and d yz magnetic moments in LnFeAsO is itiner- 
ant versus local-moment magnetisms. In addition, when 
m-pe is increased to 1.0 /ib and higher, the Fe magnetic 
moment evolves gradually to an almost spherical shape 
with no specific orbital characters. Thus, the orbital- 
distinctive magnetic moments and their effects will ap- 
pear with different strengths in various iron pnictides and 
chalcogenides for the wide range of reported niF e - 

In conventional SC, electron-phonon interactions de- 
pend on phonon polarizations, and too large interactions 
may induce structural changes that suppress SC. Like- 
wise, the Fe magnetic moment in LnFeAsO may have 
different roles in SC and AFM depending on its orbital 
character. Since our results show that the d yz charac- 
ter couples more strongly to electrons than the d xy one, 
a possibility is that the former drives AFM while the 
latter mediates SC. Another possibility is that only the 
d xy character is significant in both AFM and SC, which 
could be relevant if the d yz character has much weaker 
coupling to electrons in real samples than our results due 
to, e.g., fluctuations. An opposite scenario is that only 
the d yz character is significant in AFM and SC, which is 
less likely according to our result that SmFeAsO, having 
much higher T c than LaFeAsO when doped, shows the 
strong d xy character while LaFeAsO shows the d yz one. 

In conclusion, we have found that Fe magnetic moment 
in LnFeAsO (Ln = La to Gd) has distinctive orbital char- 
acter, d xy or dy Z , depending on Ln and uiF e . With orbital 



analysis of electronic structures, we have shown that the 
origins of the d xy and d yz orbital characters of the Fe 
magnetic moment are FS nesting and local interactions, 
respectively. These orbital characters have different cou- 
pling strength to FS electrons, so they can be identified 
by orbital-sensitive measurement of electronic structures 
and they may have different roles in SC and AFM. These 
results pave the way for understanding the role of spin 
or magnetic interactions in FeAs-based superconductors. 
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